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Abstract We report on the discovery of the first BL Herculis star displaying period- 
doubling behaviour as predicted by the hydrodynamic models of Buchler & Moska- 
lik H]. The star, with Pq ~2Ad, is located in the Galactic bulge and was discovered 
with OGLE-III photometry. We present new nonlinear convective models, which, 
together with recent evolutionary tracks, put constraints on the stellar parameters. 
In particular, we estimate the mass and metallicity of the object. 



1 Introduction 

BL Herculis stars are a subgroup of type II (or Population II) Cepheids, pulsat- 
ing with periods between 1 and 4 days (see fS] for a recent review). They show 
a singly-periodic large-amplitude light variation. Such behaviour was qualitatively 
reproduced with early nonlinear radiative models, e.g. Buchler & Moskalik HIO. 
In several of their models with periods between 2 and 2.6 days, Buchler & Moska- 
lik m found period-doubling behaviour - oscillations with periodic alternations of 
deep and shallow minima - a phenomenon not observed in any BL Her star at that 
time. Here we report on the discovery of the first BL Her star clearly showing period- 
doubling behaviour. We summarize the observations and new pulsation models for 
this star. A more detailed analysis, as well as discussion of another BL Her star in 
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Fig. 1 Pre-whitening se- 
quence forBLG184.7 133264. 
Upper panel: power spectrum 
of the original data. The 
spectrum is dominated by a 
pulsation frequency /o and its 
daily aliases. Middle panel: 
power spectrum after remov- 
ing /o and its hamionics. 
The highest peak corresponds 
to a subharmonic frequency 
j/o. The daily aliases are 
also prominent. Lower panel: 
power spectrum after remov- 
ing /o, its harmonics and its 
subharmonics. 
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Fig. 2 Light curve of 
BLG184.7 133264 phased 
with twice the pulsation pe- 
riod, 2Po. 
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which period-doubling behaviour is strongly suspected can be found in Smolec et 

al. m. 



2 Observations 



The star of interest, BLG184.7 133264 is in the Galactic bulge and was discovered 
in data collected during the third phase of the Optical Gravitational Lensing Experi- 
ment 17] . It pulsates in the fundamental mode with a period of fo = 2.4 d. The /-band 
data for the star were analysed using standard consecutive pre-whitening technique. 
Results are presented in Fig.[T] After removing the fundamental mode frequency, /o, 
and its harmonics, additional signals are visible (middle panel of Fig. [T]i. The domi- 
nant peak is located at j/o, which is a subharmonic of the primary frequency. Other 
subharmonics (j/o, 5/0 and j/o) are also present. The presence of these frequen- 
cies in the power spectrum is a characteristic signature of period doubling. In the 
time domain it means that the light curve repeats itself after two pulsation periods, 
instead of one. This gives rise to strictly periodic alternations (Fig.|2]i. 



Period doubling in BL Herculis type 



star 



3 




Fig. 3 The theoretical BL Her instability strip in the HR diagram for M = O.5OM0 and Z = 0.01 
(left panel). Nonlinear models are computed along a line of constant period, Pq = 2.4 d. The thick 
segment along this line indicates the period-doubling domain. Light curves for four selected models 
(marked with crosses) are displayed in the right panel (arbitrarily shifted in the vertical direction). 
These are to be compared with the observed light curve (represented by a Fourier fit to the data) 
plotted at the top. 



3 Theoretical models 

Using our nonlinear convective pulsation codes |I5|, we have computed several 
model sequences in order to model the period doubling in BLG184.7 133264 and to 
constrain its parameters. Computational details are presented in ||6l. Here we focus 
on presenting our best-fitting models and the resulting constraints on the parameters 
of BLG184.7 133264. 

We have considered a grid of model masses (M = O.SOMq, O.SSM©, O.6OM0 
and O.65M0) and model metalUcities (Z = 0.01, 0.001 and 0.0001). For each (M, 
Z) combination, an extensive grid of linear models was computed which covers 
the full BL Her instability strip in the HR diagram (Fig. |3]l- Next, a sequence of 
nonlinear models was computed along a line of constant period, Pq = 2.4d (right 
panel of Fig. O. Period-doubling behaviour was found over a range of luminosities. 
In Fig. [3] these are shown by the thick line for a particular (M, Z) combination. 
This domain correlates with the loci of the 3:2 half-integer resonance between the 
fundamental and first overtone modes, shown by the dashed line in Fig. [3] As shown 
by Moskalik & Buchler @, half-integer resonances are responsible for the period- 
doubling behaviour Our detailed analysis Q confirms that the 3:2 resonance indeed 
causes the period doubling behaviour in the BL Her models, as Buchler & Moskalik 
have already shown 

The model light curves were compared with the observations. Our best-fitting 
models closely match the pulsation amplitude and also match the amplitude of the 
alternations (i.e., the amplitude of the subharmonic peak in the frequency spectrum) 
of BLG184.7 133264. A good match is possible only for models with the highest 
metallicity, Z = 0.01. For lower metallicities, the pulsation amplitudes are lower 
than observed and, in addition, the ampUtudes of the alternations are much larger 
than observed. This result is independent of the adopted values of the convective 
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Fig. 4 The BaSTI (H hor- 
izontal branch evolutionary 
tracks for Z = 0.01. Evo- 
lutionary tracks start at the 
Zero-Age Horizontal Branch 
(solid, horizontally running 
line). Each track is labeled 
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with the shaded area indicat- 
ing the BL Her domain with 
fundamental mode periods 
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parameters, as analysed in detail in 161. The three best-fitting models, all with Z = 
0.01, have different masses: M = O.SOMy , O.SSM.j and O.6OM0. Based on pulsation 
computations alone, we cannot decide on the best model. However, since all these 
models fall roughly in the same place in the HR diagram (see Fig.|4|i, evolutionary 
tracks may provide further constraints. 

In Fig.|4]we plot horizontal branch evolutionary tracks from the BaSTI database 
corresponding to Z = O.OI ||4|. It is clear that only our least massive model (M = 
O.5OM0) fits the evolutionary scenario. The tracks for larger masses run well beyond 
the instability strip. The light curve of the best model is shown in Fig. [3] (second 
model light curve from top, plotted with a thick line). 
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